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In most cases, during the unsteady state, more of the heat arriving 
at the surface of the liquid droplet goes to sensible heat than goes to 
latent heat of vanorization. This is esnecially true for low-volatility 
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mass transfer decreases with time, although the droplet is being heated 
up. The time required by a certain size droplet to reach its wet-bulb 
temperature is less for higher volatility fuels at the same ambient and 
initial conditions . 

It should be realized that the calculated curves depend upon the 
use of heat and mass transfer correlations that are extended far beyond 
the conditions under which they were obtained. Thus, extensive use of 
the values shown on the calculated curves should not be made until either 
the curves or the extrapolations of the heat and mass transfer correla- 
tions have been checked experimentally. 

In addition to the theoretical work, in order to verify the results 
of the calculations, which showed the unsteady state to be of importance, 
apparatus was set up to obtain experimental temperature-time histories 
of comparatively large -Iroplets vaporizing in air streams at constant 
velocities and atmospheric pressui'e. This apparatus is described and a 
few experimental e a^ves are included. These curves show the relative 
distribution of the time Interval between the uj'isteady and steady states 
when different air temperatures and different fuels are used. 


INTRODUCTION 


Knowledge of the phases through which a heterogeneous mixture of 
air and fuel passes before ignition taJces place is essential for the 
rational design of jet-engine combustion chambers. For discussion pur- 
poses, the time elapsed before ignition is often divided into physical 
and chemical components, although the relative magnitudes of these two 
components are not known. Since the physical component inevitably pre- 
cedes and perhaps in part overlaps the chemical component, knowledge of 
the magnitude of and the factors affecting trie ^ lysical component is of 
obvious importance in any study of the time inturva^l before ignition. 

Fuel injected through a nozzle into a combustion chamber is believed 
to pass through the following stages: The fuel leaves the nozzle orifices 

as a ligament or sheet. This ligament or sheet then breaks down into dif- 
ferent size droplets. Appreciable heat transfer from the air present in 
the combustion chamber to the fuel occurs only after the ligament has 
broken up into droplets. This is due to the appreciable increase of sur- 
face area of the fuel that is exposed to the air. As heat transfer takes 
place, the droplets heat up (or possibly cool dov, . depending upon ambient 
air conditions and initial fuel tempe -’ature) and a the same time lose 
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relative to the air by aerodynamic drag forces. After a certain time 
has elapsed, each droplet attains an asymptotic or equilibrium tempera- 
ture equal to the wet-bulb temperature corresponding to the conditions 
present at that moment . 

The larger droplets are slower in attaining equilibrium conditions; 
but, although they have essentially the same initial speed as the smaller 
droplets, they are slowed down relative to the air at a lesser rate, and 
they move along with the air ahead of the smaller ones that were injected 
at the same instant. The smaller drops, however, give away their mass 
and completely vaporize faster and travel a shorter distance through the 
combustion chamber than the larger droplets. 


A cloud of vapor due to the smaller droplets is thus rapidly formed 
and moves along with the air . The mass of vapor given away by the 
incoming larger droplets is added to this vapor cloud. Somewhere in the 
combustion chamber a combustible mixture of air and fuel vapor is formed 
and is ready for ignition at that point. In the absence of outside 
ignition, this mixture must be heated by the air until it reaches its self- 
ignition temperature if combustion is to take place. By definition, the 
physical portion of the ignition delay extends until a combustible mixt\xre 
is formed and heated to the ignition temperature. Since some chemical 
reaction may take place at temperatures lower than the ignition tempera- 
tirre, the chemical portion of the ignition delay period probably overlaps 
the physical portion. The vaporization process is thus an essential part 
of the ignition delay period. 

Although the fuel droplets reach their equilibrium or wet -bulb tem- 
peratures asymptotically with time, the vaporization process can roughly 
be separated into the unsteady state and the steady or equilibrium state. 
Numerous investigations have been conducted covering the latter stage 
(refs. 1 and 2), but only a few attempts have been made to study the 
imsteady state (ref. 5). It was the intent of this investigation to 
determine theoretically temperature-time histories of droplets of dif- 
ferent size and fuel composition under different ambient air conditions 
in an attempt to determine whether the lansteady-state period is of 
importance in the history of the fuel droplets previous to ignition. 


This report, then, presents a detailed study of the unsteady-state 
portion of the total vaporization time of single fuel droplets injected 
into air. Calculations were performed for numerous combinations of the 
following conditions: 


Fuels 

Initial radii, p 

Initial droplet temperatures, 9F 
Initial droplet velocities, fps 

Air pressures, atm 

Air temperatures , °F 


n-octane, benzene, n-pentane 

10, 50, 100 

50, 200 

100, Aoo 

0 . 5 , 1 , 5 

100, 500, 1,000 
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The three fuels have standard atmospheric boiling temperatures of 96.9° 
176.2°, and 258.2° F, respectively. The calculations covered all com- 
binations of the above conditions except when the initial droplet tem- 
peratures were higher than the boiling points of the fuels at the corre 
spending air pressures . 

The properties of the three fuels under investigation and for air 
were obtained from various references as shown in table I. 

This work was conducted at the University of Wisconsin under the 
sponsorship and with the financial assistance of the National Advisory 
Committee for Aeronautics. 


SYMBOLS 

Ac cross-sectional area of liquid droplet, sq in. 

Aq surface area of liquid droplet, sq in. 

Bo thickness of air-vapor film sxirrounding droplet, in. 

Cp coefficient of drag for spheres, unitless 

specific heat of air at constant pressure, Btu/(lb)(°F) 

Cp-p specific heat of fuel vapor at constant pressure, 

Btu/(lb)(°F) 

Cp^ specific heat of liquid fuel, Btu/(lb)(°F) 

Cpjj^ average specific heat of air -vapor mixture in film at 

constant pressure, Btu/(lb)(^) 

Bv diffusion coefficient of air-vapor system, sq in. /sec 

d penetration or distance traveled relative to air, in. 

F drag force, lb 

Sa weight fraction of air in film, unitless 

8 f weight fraction of fuel vapor in film, unitless 

h coefficient of heat transfer in film in absence of mass 

transfer, based on surface area of liquid droplet Aq, 

Btu/(sq in. ) (sec) (°F) 



NACA TN 5179 


5 


K^, K3, 

ka 

V 

M 

Ma 

M-p 


H 


'■m 


m 


N 


Nu 


N- 


Nu 


N- 


Pr 




Re 


N, 


Sc 


a 


PaB 


coefficient of heat transfer to liquid^ Btu/(sq in,)(sec)(^) 
. . . constants 

thermal conductivity of air^ Btu/( in. ) ( sec) (^F) 

thermal conductivity of fuel vapor ^ Btu/(in. ) (sec) (^F) 

coefficient of mass transfer in film^ l/sec 

thermal conductivity of liquid fuel_, Btu/(in. ) (sec) (^F) 

average thermal conductivity of air -vapor mixture in film^ 
Btu/( in. ) (sec) (^F) 

ratio of instantaneous mass to original mass of droplet^ 
percent 

molecular weight of air^ lb /mole 

molecular weight of fuel vapor^ Ib/mole 

apparent molecular weight of air-vapor mixture^ lb /mole 

mass of liquid droplet, lb 

Nusselt number for heat transfer, unitless 

Nusselt number for mass transfer, unitless 

Prandtl number, Cp unitless 

Reynolds number, 2rQU%^nm> unitless 
Schmidt number Pi^i/^m^v^ unities s 
total pressure, Ib/sq in. abs 

average partial pressure of air in film, P^i - 
Ib/sq in, abs 

partial pressure of air at outer edge of film, 

Ib/sq in. abs 
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partial pressure of air at liquid droplet surface, 
lb/ sq in. abs 

average partial pressure of fuel vapor in film, 

Ib/sq in. abs 

partial pressure of fuel vapor at outer edge of film, 
considered to be zero, Ib/sq in. abs 

partial pressure of fuel vapor at liquid droplet surface, 
equal to vapor pressiire of fuel at droplet ten5>era- 
ture T^, Ib/sq in. abs 

total heat transfer from air to droplet, Btu/sec 

sensible heat received by droplet, Btu/sec 

heat carried away from droplet by diffusing vapor in form 
of superheat, Btu/sec 

heat received at droplet surface, Btu/sec 

heat of vaporization, Btu/sec 

universal gas constant, in-lb/(mole) (°F) 

radius at any point in film, in. 

molal rate of diffusion of air in equimolal diffusion, 
moles/(sec) (sq in.) 

molal rate of diffusion of fuel vapor in equimolal 

diffusion, based on surface area of liquid droplet Aq, 
moles/(sec) (sq in.) 

radius of liquid droplet, in. 

temperature in film at raditis r, °R 

asymptotic or wet-bulb temperature of droplet, °R 

temperature of air outside film, °R 

temperature of liquid droplet, 

logarithmic mean tengjerature in film, °R 
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u 

w 

X 


y 


z 


a 


0 

A 


Pa 

PL 


droplet velocity relative to air, in* /sec 
mass of fuel vapor diffused out, Ib/sec 
proportionality factor 


fraction of film thickness at radius r, 



unitless 


factor, unitless 

correction factor of mass transfer due to semipermeahility 
of film, unitless 

time, sec 

latent heat of vaporization of fuel, Btu/lh 
viscosity of air, lb/(in.) (sec) 

viscosity of fuel vapor, lb/(in. ) (sec) 

average viscosity of air- vapor mixture in filin, lb/( in. ) ( sec) 

dens i ty of a i r ^ lb / cu in . 

density of liquid droplet, lb/ cu in. 

average density of air-vapor mixture in film, Ib/cu in. 


ANALYSIS 


In the following theoretical study of the vaporization of single 
droplets, two assumption:, were made. The first assumption was that an 
endless supply of air Siu-rounded the droplet at all times. This means 
that there is no interaction between different droplets in the same air 
stream, or that only a single droplet was considered present. The 
partial pressure due to i.he fuel vapor at the i. t.er edge of the air- 
vapor ' .^urroi.. : ia< the m , et was thus considered to be zero at all 
times. Tlie seconu H.ssumpriAjn was that the air vessiure and temperature 
rertiained nnc narigovi t.r.r,- ■■.-■-t the vaporization process. It kiiowu 



8 


NACA TW 5179 


that during injection in a continuous -flow air stream the air temperature 
decreases because of heat transfer and that in the absence of friction 
its pressure should increase because of the slowing down of the stream; 
the second assumption is consistent, however, with the assumption of a 
single droplet in a large volume of air. 

Figure 1 shows a spherically symmetrical droplet surrounded by a 
film containing an air - fuel-vapor mixture. At any one instant, the 
mixture strength varies from a maximum at the droplet surface to a mini- 
mum at the outer edge of the film. The shape of the temperature gradient 
curve depends on the relative temperatures of the fuel droplet and air 
and the mass of vapor being diffused out at that instant. Under the 
ambient conditions considered in this study, eliminating the case where 
the droplet was introduced with a vapor pressure higher than the total 
pressure, the vaporization process can take place in three different 
cases: 

(1) When the initial temperature of the droplet is lower than the 
air temperature and lower than its own wet-bulb temperature, in which 
case the droplet heats up continuously during the unsteady-state portion 
of the vaporization time. This case cort^jrises the larger part of the 
calculated conditions contained in this report. 

(2) When the initial ten5)erature of the droplet is lower than the 
air temperature but higher than its own wet-bulb temperature at the 
ambient air conditions. In this case the droplet cools down continu- 
ously while vaporizing and the sensible heat lost by the liquid supplies 
part of the heat of vaporization of the diffusing vapor. Heat is always 
transferred from the air to the droplet, however. 

(3) When the initial temperature of the droplet is higher than both 
the air temperature and its own wet-bulb temperature. In this case the 
droplet cools down continuously until it reaches its wet-bulb temperature 
which is lower than the air temperature. During that time when the droplet 
temperature is higher than the air temperature, there is a net heat trans- 
fer from the droplet to the air and thus heat and mass transfer take place 
in the same direction. When the droplet temperature has dropped below 

the air temperature, the process is the same as that in case (2). 

It is convenient to divide the analysis into three parts - the heat 
transfer, the mass transfer, and the velocity change. 


Heat Transfer 

Referring to figure 1 , and assuming that the droplet is vaporizing 
under the conditions of case (l), the total or net heat transfer Q from 
the air to the film surrounding the droplet, prior to ignition, goes three 
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ways: (l) To heat up the liquid droplet, ( 2 ) to vaporize the 

liquid, Q.; and (5) to be carried or peeled bach with the diffusing 
vapor in the form of superheat, Qs • The heat that arrives at the drop 
let surface is denoted by and equals the sum of and Q^. 


It should be noted that In developing the relationship used Q does 
not include radiant heat transfer. 


If at any 
at a radius r 
considered: 


instant during vaporization a point in the air-vapor film 
from the center of the droplet and at temperature T is 


^ = Ql + 


= Q - Qg 


( 1 ) 




= hA, 


o f 


( 2 ) 


where w 


is the rate at which the vapor is diffusing out, ^ 
xs bne xau ... - h is the coefficient of 

is 


fraction of the film thickness at radius r^, 

heat transfer through the film in the absence of mass transfer, Cp^ 

the specific heat of the fuel vapor, T is the terr^jerat^e of the mix- 
ture at the position y, and Tl is the temperature of the surface of 

the droplet. 

Rearranging equation (2) gives: 


dy = 


hA. 




dT 


'Pf 


Integrating across the film between y =0 to y= 1.0 where T 
to T = Tg and considering the specific heat of the fuel vapor Cp^ 

to be constant throughout the film: 


hA. 


1 = 


log. 


+ “Cpf(TB - Tl) 




Civ 
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or by rearranging and expressing as an exponential: 

^wCp^^hAo _ - Tl) 

% 


Defining z 


as wCp^^hAo 


gives: 





a, = m„(Tj - 

- 1 


(3) 





(^) 


Thus the factor z /(e^ - l) represents that fraction of the total 
heat transfer Q from the air that finally arrives at the surface of 
the liquid droplet and supplies the latent heat of vaporization for the 
vapor diffusing out as well as the sensible heat added to the liquid 

droplet itself. The factor zj{e^ - l) therefore represents a correc- 
tion factor to the heat transfer coefficient h without mass transfer 
from the droplet. 


Several investigators have presented correlations for the heat trans- 
fer coefficient h. However^ as will appear later, it was necessary for 
these studies that a mass transfer coefficient also be available. The 
correlations used in this paper were taken from Ranz and Marshall's work 
(ref. 1). It was felt that the nature of the investigation conducted in 
their work was similar to the problems being studied in this paper, 
although the maximum air temperatures used in their work were not so high 
as the maximum air temperatures used in this investigation. 

From their correlation the Nusselt number for heat transfer used is: 


W- 


Nu = 


h(2ro) 
\ 


2 + 0.6(Np^)V5jrjj^^ji/2 


(5) 


NACA TN 5179 


11 


The factors in this equation were evaluated for the average film con- 
ditions. This required knowledge of the air and droplet temperatures 
as well as the partial pressures in the film (see appendixes A and B) . 
The Ran?, and Marshall correlations were obtained from experiments con- 
ducted at low temperatures and low vapor pressures. In their work, 

therefore, the correction factor z/(e 2 - 1 ) was unity and it is not 
clear whether their correlation includes this factor as a variable. 

Since the correlation was to be used under conditions where this factor 
was not imity, and since the correlation is similar in form to those 
obtained without mass transfer, it was decided to extrapolate the corre- 
lation by assuming that it did not include the correction factor. The 
correction factor therefore was computed and used in all of the calcu- 
lations presented herein. 

It can be seen from the above equations that a calculation of the 
heat transfer through the film at any instant during the vaporization 
process of the droplet would require the knowledge of the temperatures 
on both sides of the air-vapor film sin-rounding the droplet, of the 
velocity of the droplet relative to the air, and of its radius, which 
is dependent on the mass transferred up to that point and on the droplet 
temperature . 


Mass Transfer 

Neglecting thermal diffusion and assuming all diffusion to result 
from the driving force of a concentration or partial pressinre gradient 
which exists in the direction of diffusion, the following equations are 
obtained for counterdiffusion systems (ref. 4 ): 


dpj. 

dr 


- '■aPf) 


( 6 ) 


and 


_ (rt)^ 

xPrp 


(7) 


In this case 
total pressure is 


equimolal rates of diffusion must 
to he maintained, or: 


take place if the 


J'f = -^a 


( 8 ) 
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and 


= Pa + Pf (9) 

However, in most cases of a fuel droplet vaporizing in air essen- 
tially unidirectional diffusion takes place. Equimolal diffusion would 
be approached only if a droplet approached its critical temperature 
during the heating-up period. In the cases under consideration in this 
report essentially unidirectional diffusion exists; that is, there will 
be diffusion of the fuel vapor away from the droplet while the rate of 
diffusion of the air in the opposite direction, toward the droplet, is 
zero. In this case. 


^a = 0 

Combining equations ( 6 ), (7)^ and (lO) and rearranging; 


(10) 


_ RT 
dr D^,j, 


( 11 ) 


and 


Integrating between r = 
Pa = PaB‘- 


and r = r^ + Bq where Pg 


PaL 


rf = 


RTB^ 


log^ 


^aB 

PaL 


Dy 

RTBq 




T 


PaB " PaL 


log. 


aB 


■aL 


(12) 


This is the equation for the rate of diffusion of the fuel vapor in terms 
of the partial pressures of the air. Using equation (9) and assuming that 
the partial pressure of the fuel vapor at the outer edge of the film pj^ 

is zero, an assumption discussed earlier in the analysis, the rate of dif- 
fusion of the fuel vapor can be written in terms of partial pressures in 
the film as follows: 
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= 


RTB, 


PfL 


o 


PfL 


log, 


■T 


P - P 
T ^fL 


(13) 


where 


= 


Dy 

RTB,- 




a 


PfL 


log. 


T 


Pp - PfL 


(li+) 


Equation (l 4 ) states that the unidirectional diffusion of the fuel 
vapor into the air equals the equimolal rate of diffusion^ which is a 
function of the concentration, times the factor a which is greater than 
unity. This increase in rate represents a correction for the fact that 
the net rate of movement of the air is zero, while the diffusion veloci- 
ties of the fuel vapor and the air relative to each other must he 
maintained . 

Because of the difficulty in evaluating Bq and by analogy with 
heat transfer a semiempirical, semitheoretlcal relationship is set up 
for mass transfer. Thus it is customarily considered on a semiempirical 
basis that the mass of fuel vapor diffused out per unit time under the 
conditions of interest here is given by 

^ = VgPfL°- 

where kg is the mass transfer coefficient for the film in the case of 

equimolal diffusion. Ranz and Marshall (ref. l) present a correlation 
for kg, but again it was necessary to extrapolate their correlation to 
conditions far removed from the experimental conditions imder which it 
was obtained. In making this extrapolation the line of reasoning was 
as follows: The correlation was for the film. Therefore, when the 

density of the air and of the film were markedly different, the mean 
density of the film should be used instead of Another choice 

to be made in the extrapolation was whether the Ranz and Marshall corre- 
lation is for kg or for the product of k and a, since the data for 
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the Ranz and Marshall correlation were taken under conditions where a 
was essentially unity and therefore not a variable. The factor p is 

sometimes used as a correction factor instead of a (ref. 4 ). However 
when it is so used^ it is a logarithmic mean average, while in all of 
their calculations for their correlations Ranz and Marshall used an 
arithmetic mean average. The difference between the two averages is 
almost undetectable over the range of conditions covered in their experi- 
ments and thus it is not clear whether their correlation included a 
via the factor It is also not clear whether the use of the arith- 

metic mean rather than the logarithmic mean is an essential part of their 
correlation. For extrapolation purposes for these computations, it was 
decided to use p^^ as an arithmetic average and to assume that a was 

not included in the correlation of Ranz and Marshall. Thus the final 
equation used was 


%u’ 


g^oPa^g 

®vPm 


2 + 0.6(Ng^)V5(KR^^l/2 


(16) 


factors in the mass transfer equations were also computed for 
the average film conditions (see appendix B). The values of the dif- 
fusion coefficient Dy were arrived at by using the techniques recom- 
mended by Hirschfelder, Curtiss, and Bird (ref. 5). It can be seen again 
from the above equation that, to calculate the mass transfer at any 
instant, a knowledge of the heat transfer and of the radius and velocity 
of the droplet at that instant is essential. 


Velocity Equation 

The velocity of the droplet at any instant during the vaporization 
process is essential to the calculation of both the heat and mass trans- 
fer at that instant, since it enters into both correlations. The droplet 
is assumed to be injected with a certain initial velocity relative to 
the air. Aerodynamic drag forces will then either slow it down or speed 
it up so that its velocity approaches that of the air. The less the 
density of the air, that is, the lowei' its pressm'e and the higher its 
temperature, the less effective the drag forces and the greater and more 
rapid the penetration through the combustion chamber. 

From aerodynamic theory the drag force exerted by a moving fluid 
upon an immersed body may be calculated from the equation 


= CdAcP^ ^ 


(17) 
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where Cj) is the coefficient of drag for the given body form and is a 
function of the Reynolds number; the cross-sectional area of the 

body; and the density of the fluid. 

For a spherical droplet: 


or 






^ - 1 Cj) — — 

de 8 ® Pl 


( 18 ) 


Here Cj) is taken as the drag coefficient for spheres (ref. 6) at 
the Reynolds number for which the calculation is performed. It has been 
pointed out by Hughes and Gilliland (ref. 7) that the drag coefficients 
for solid spheres and for droplets are not equal. However, the condi- 
tions for which the corrections are shown do not correspond to the con- 
ditions for which the calculations were performed and therefore no 
correction was made. Firrthermore, extrapolation of the data of Hioghes 
and Gilliland would indicate that the correction should be close to 
unity for the conditions considered here. It is again evident that, to 
calculate the deceleration of the droplet at any instant, a knowledge 
of both the heat and mass transfers is essential. 


Method of Calculation 

The dependency of the heat transfer, mass transfer, and velocity 
equations upon each other and the complexity of the equations neces- 
sitated a stepwise integration process where the values required are 
evaluated for small increments of droplet temperature over a small 
period of time during which all properties of the droplet and fluid may 
be assumed to be constant. A sample calculation is shown in appendix A. 
The procediire followed was first to assume an increment of time A0, a 
final velocity U 2 ^ and a final radius r ^2 for a certain increment of 

droplet temperature the initial values of velocity Up, radius ^q1> 

and temperature T^p for that increment being already known from the 
previous step, or, if this were the first increment, being specified. 

For that increment the mean values of the droplet temperature T^, 
radius r^, and velocity U are calculated. Knowing the pressure and 

temperature of the air, which are assumed to remain constant throughout 
the process, the various properties in the film are calculated (see 
appendix B) . 
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Then 


is calculated from equation (l 8 ) and U2 is evaluated: 


U2 = A6 (19) 

This is repeated several times, correcting for each time, until 

a value of U2 is found that satisfies the velocity equation and the 
Reynolds number is fixed at the mean velocity of the increment U. 

The next step would be to calculate the mass transfer coefficient kg. 

o 

and the average mass transfer w in pounds per second for the increment 
from equations (15) and (16). This requires a knowledge of the partial 
pressure of the fuel. The total mass transfer for the increment is then 
found from w AB and the average droplet mass m would be 


^1 

m = 

2 


(*^1 ” ^ 
2 


(20) 


In this equation m^^ would be specified or known from calculations of 
the previous step. 

The heat transfer to the droplet surface is then calculated 

from equations (3) and ( 4 -) and the remaining calculations take place as 
follows : 


Ql = Qv - 

For all calculations shown here it was assumed that the thermal conduc- 
tivity of the fuel was high and the droplet temperature was uniform. The 
validity of these assumptions is discussed in appendix C. Therefore, 


A 0 = 


mCp ATl 


Qt 


3m2 


■o2 






( 21 ) 


1/3 


(22) 
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The values of U 2 , A0, and calculated from eqmtions ( 19 ), 

( 21 ), and ( 22 ) are compared with those assumed at the beginning of the 
computations and are usually recalculated at least once more. The final 
values for this increment then become the initial values of the next 
increment and the computations are continued for small increments of 
ZYTl until Ql is found to approach zero. At this point the temperature 

asymptote is about to be reached. Here the computations were stopped, 
since only the unsteady state was being studied. Since the asymptotic 
temperature was not known for most conditions prior to each calcula- 
tion, /:^rL and the number of steps were difficult to determine. It was 
found that between 12 and I 6 steps usually gave sufficient accuracy. 

A model 1, IBM Card-Programmed Electronic Calculator was used for 
the computations . The special program to solve this problem uses the 
floating-decimal general-purpose connections of the calculator. One 
complete calculation for specified droplet conditions of initial radius, 
temperature, and velocity and air conditions of pressure and temperature 
takes approximately 2.0 hours on the CPC not counting the setup time, 
computation of the properties at the successive droplet temperature incre- 
ments, and card prepunching. 

A typical vaporization history calculated by this method is shown 
in figure 2 for a droplet of n-octane of 50 -micron initial radius, 

50 ® F initial temperature, and 100-fps initial velocity with respect to 
air, vaporizing under air ambient conditions of 0 . 5-a-'tmosphere pressure 
and 1,000° F temperature. 


Steady State 

The droplet temperature approaches the wet-bulb temperature asymp- 
totically dixring the vaporization process. At the wet-bulb or asymptotic 
temperature 

Qv = (23) 

Substituting values from equation (j) and changing to there 

is obtained: 




( 24 ) 
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Substituting wCp^yiiA^ for z and rearranging: 

^ (e" - 1) (25) 

Cpf 

This is a simple equation for the asymptotic or wet-bulb temperature 
as a function of the air temperature and the mass and heat transfers, A 
knowledge of the latter two values is essential for a solution of this 
equation. It can^ however, be rewritten using only the properties of 
the liquid and film as follows: Combining equations ( 24 ) and (15) and 

rearranging give: 




as 


PfpA 


a 


h 


- 1 


Combining the above equation with equations (5) and (16) and rearranging: 


T, 


as 


= Tx 


PyPmPfL^ 

PaV 



(26) 


Solution of equation (26) must necessarily be by trial and error^ 
since the component properties are all functions of both and Tn 

as well as of the air pressure and the fuel used. The factor zj(e^ - l) 
must be approximated since this factor is affected by the volatility of 
the fuel. It can be stated that for very low volatility fuels this factor 
can usually be taken as luiity with little error. For high- volatility 
fuels at low air temperatures and high air pressures z/(e^ - 1) can 
still be taken as unity with little error. For example, using n-pentane, 
this factor has a value of 0.957 s-t 100 *^ F air temperature and 5 -atmosphere 
air pressiure while n-octane, a less volatile fuel, has a value of O.998 
for this factor at the same conditions. However, for high-volatility 
fuels at high air temperatures and low air pressures lower values have to 
be used. n-Pentane has a value of about O.65 at 500 ° F and 1 atmosphere 
air conditions . 


The mass transfer during the steady-state period of vaporization can 
be determined by combining equations (15) and (16) as follows: 


2D^PmPfL^ 
w = r 

P.. 


o 


2 + 0.6 
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Since the temperatures across the film have become essentially constant^ 
the above equation can be rewritten as: 


w = K,r^[i + K2(r^u)l/g 


(27) 


where Kq and K2 are constants . 


Equation (27) shows that for a particular fuel the rate of nkuss 
transfer during the steady state is a fimction of only the radius and 
velocity of the droplet at any instant, provided the air temperature 
and pressure remain unchanged. 


Equation (26) for the steady state where the temperatures across 
the film are fixed, and therefore all the properties in it, can be 
rewritten as 


as 


= T- 


B 


K, 


2 t 

2 + K5(roU)l/2 


(28) 


where and are constants for any particular fuel and air 

conditions at the asymptotic or wet-hulh temperature Therefore^ 

depends upon the ratio of the Nusselt numbers for mass and heat 

transfer. It can he deduced^ after examining the above equation^ tliat 
the wet-bulb temperature is dependent on the droplet radius and velocity. 
At high droplet speeds or for relatively large droplet radii the fac- 
tor 2 in both the Nusselt numbers becomes relatively negligible and the 
ratio of the Nusselt numbers and consequently the wet-bulb temperature 
become essentially constant. However, at lower values of droplet speed 
or droplet radius the wet-bulb temperature rises slightly until either 
the speeds or radii approach zero, where the Nusselt number ratio 
approaches unity and the wet-bulb temperature is again fixed. 

This dependence proved to be small for most conditions of interest 
in this paper. Figure ^(si) shows ATg^g as calculated for various veloci 

ties, where ATas defined as the deviation of the actual wet-bulb tern 
perature from that calculated at large Reynolds numbers; figure 5 (b) 
shows ATas for various radii for n-pentane at 1-atmosphere air pressure 
and 100 ^ F air temperature. It should be noted that no attempt was made 
to include a dependence of the vapor pressure of fuels on dropleh size 
as the droplet radius became very small in calculating the information 
shown in figure 3(b) . It is to be expected that the rise of vapor pres- 
sure as corrected for very small droplet radii would lower the values of 
the wet-bulb temperature again at small drop radii. 
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Further examination of equation (26) reveals that the wet-bulb tem- 
perati^e approaches the air temperature in case of low-volatility fuels. 
This is to be expected, since the mass transfer has to be high enoiigh to 
carry away the entire heat transferred to the droplet. In the case of 
no mass transfer, the wet-bulb temperature equals the air temperature. 


APPAEATUS 


After investigating theoretically the importance of the imsteady 
state during the vaporization period and conducting a number of calcu- 
lations in the unsteady state, it was felt that experimental verifica- 
tion of its importance was desirable. Accordingly, experimental appara- 
tus has been set up in an attempt to establish the Importance of the 
unsteady state and a few experimental records have been obtained. 

Figure 4 shows a schematic diagram of the apparatus constructed to 
determine the temperature-time histories of vaporizing droplets. Air 
is tapped from the house line at 100 Ib/sq in. abs and passed through 
a porous-stone filter to remove water, oil, and other Impurities. Its 
pressiire is then controlled by a pressure regulator, after which it 

passes through a 2- -inch pipe where its pressure and temperature are 

measured by a precision Bourdon gage and a thermometer. The diameter of 
the pipe at this measuring section is made large enough so that the dif- 
ference between total and static temperatures is negligible at the highest 
rate of flow to be used. The air is then passed through either the 
5/32- or the 5/l6-inch-diameter critical-flow orifice. The air is then 
heated with an electrical resistance heater, after which it passes through 
a calming section to a l4-inch-long section of 2-inch flexible tubing 
which terminates in a 1-inch I.S.A. standard flow nozzle. The exit sec- 
tion of the nozzle is covered with a l60-mesh screen. The screen assists 
in eqijalizing the velocity over the nozzle opening to a value calculated 
from the average volumetric velocity. 

The liquid droplet is hung on the junction of the thermocouple 
about 1/2 inch above the exit section of the nozzle. A 0.010- inch iron- 
constantan thermocouple was used for the preliminary investigations 
reported herein. A clean and small thermocouple Junction was found essen- 
tial for a minimum of droplet distortion. It was found that droplets of 
TOO- to 900-micron radius were the best size to be used with the above 
thermocouple wire size without distortion during the major portion of the 
heating-up and vaporization period. 


An illuminating system was provided for viewing the droplet. A 
500-watt General Electric TSi projection lamp was used. A heat-absorbing 
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glass was placed in the path of the light to minimize the heating of 
the droplet by radiation from the lamp. An optical system having a maxi- 
mum magnification of 86 was provided for measuring the diameter of the 
droplet. The optical system was calibrated by placing a wire of known 
diameter in place of the droplet. 

The experimental procedure was as follows: The air was heated to 

the desired temperature with the exit nozzle swung away from the 
thermocouple . 

The exit velocity of the air stream was adjusted by varying the 
inlet air pressure ahead of the critical flow orifice. Low air veloci- 
ties were used to avoid blowing the droplet off the themocouple junction. 

The droplet was then formed and suspended on the thermocouple using 
a syringe and hypodermic needle. The initial droplet radius was then 
measured by the optical magnification system. The air nozzle at the 
end of the flexible tubing section was then suddenly swung under the 
droplet. The output of the thermocouple was measured and recorded by 
a 1 -second full-scale model G Leeds & Northrup electronic self-balancing 
Speedomax potentiometer. The temperature of the droplet was measured by 
the calibrated thermocouple output, while the time scale was computed 
from the known speed of the recorder chart. 


RESULTS AM) DISCUSSION 


The typical calculated vaporization-time histories shown in fig- 
ures 2 and 5 through 8 contain the following information: 

(1) Droplet temperature, Tl, ^ 

(2) Droplet mass, percent of original, M, percent 

(3) Radius, r^, inches 

( 4 ) Velocity relative to the air, U, fps 

(5) Penetration or distance traveled relative to the air, d, inches 

(6) Instantaneous rate of mass transfer, w, pounds per second 

(7) Instantaneous total heat transfer from the air to the film, Q, 

Btu/sec 

(8) Instantaneous heat of vaporization, Q^, Btu/sec 

(9) Instantaneous heat transfer to the liquid droplet, Qj^, Btu/sec 

The fuel, initial droplet radius, initial droplet temperature, initial 
droplet speed relative to the air, air pressure, and air temperature are 
indicated in the title of each figure. The portion of the total heat trans- 
fer Q that is carried back with the diffusing vapor in the form of 
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superheat can be determined from the curves by taking the difference 
between Q and the SLun of and Q^. 

Figures 2 and 5 through 8 represent the results of only 7 out of 
the 252 eases for which calculations were performed. A complete set of 
graphs presenting the results of all the calculations is contained in 
the manuscript copy of this report. This copy is available for loan or 
reference in the Division of Research Information, National Advisory 
Committee for Aeronautics, Washington, D. C. 

Examination of these typical graphs (figs. 2 and 5 to 8) yields the 
following conclusions: 

For low-volatility fuels^ for example^ figure 2 for n-octane^ 
speaking of the heat transferred to the surface of tne liquid droplet, 
the amount of heat going to sensible heat far exceeds the amount of heat 
going to supply the latent heat of vaporization during most of the unsteady- 
state or heating-up period and especially right after injection. The 
higher the air temperature or the higher the air pressure, the greater the 
portion going to sensible heat. Only as the droplet approaches the steady 
state does the relation reverse itself, as required by the fact that 

becomes negligible and essentially equal to zero during the steady state. 
This effect is, of course, largely dependent on the initial injection 
temperature of the liquid fuel. At high initial liquid temperatures sen- 
sible heat can be lost by the droplet or a net heat transfer Q from the 
liquid to the air may even occur (e.g., fig. 5 for n-octane at 200° F 
initial liquid temperature and 100^ F air temperature) 

For high-volatility fuels, however (e.g., fig. 6(a) for n-pentane 
at 1-atmosphere air pressure and 500° F air temperature) , the latent 
heat of vaporization exceeds the amount going to sensible heat. In some 
cases, especially at low air pressures (e.g., fig. 6(b) for n-pentane), 
the liquid undergoes a cooling down or loses sensible heat at 0. 5 -at mo sphere 
air pressure regardless of the air temperature. In this case the net heat 
transfer from the air goes entirely into latent heat of vaporization. The 
balance of the heat required to supply the latent heat of vaporization is 
drawn from the droplet itself . Only at high air temperatures and high air 
pressures does the sensible-heat requirement of the droplet exceed that of 
the heat of vaporization, and only for a shorter portion of the unsteady 

The rate of mass transfer usually increases with time as the droplet 
is heating up. The increase is much more pronounced at higher air tem- 
peratures and lower air pressures (e.g., fig. 2 for n-octane at 
0.5-atmosphere air pressure and 1,000° F air temperature). At high air 
pressures and low air temperatures (e.g., figs. J and 8(a) for benzene 
and n-octane at 5-atmosphere air pressure and 100° F air temperature) 
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the rate of mass transfer undergoes a decrease, although the droplet is 
heating up. The reason for this is that in this case the denser air 
slows down the droplet very rapidly with a consequent rapid drop in the 
Reynolds number which offsets the favorable effect of the exponential 
increase in the vapor pressure of the liquid fuel with temperature on 
the rate of mass transfer w. At intermediate cases, a rise in w 
followed by a drop can be noticed (e.g., fig. 8(b) for n-octane at 
1-atmosphere air pressure and 100° F air temperature). In the case of 
droplets cooling down during the unsteady state, the rate of mass trans- 
fer consistently decreases with time. In the steady state the rate of 
mass transfer becomes a function of only the radius and velocity of the 
droplet and thus undergoes a small but steady decrease in value (see 
eq. (27)). 

For comparison purposes the three-dimensional graphs shovn in 
figs. 9 through 12 are presented. These are temperature -time history 
curves for the three fuels at different fuel and ambient air conditions. 

In all these curves the initial temperatirre of the droplets was 50° F 
and the initial velocity relative to the air was 100 fps . 

Figure 9(a) shows temperature-time histories of n-octane at an air 
temperature of 500° F, for various droplet sizes and air pressures. It 
can be seen that for the same air pressures but different radius droplets 
the asymptotic temperatixres are eq.ual and the three curves representing 
the three sizes at any one air pressure form a horizontal plane after all 
three droplets have reached their steady states. The ends of the hori- 
zontal planes shown by the broken lines have no particular meaning and 
were drawn for illustrative purposes only, since the 10-micron droplets 
would vaporize before reaching this point and the 100-micron droplets 
would last longer. Figure 9(^) shows similar temperature-time histories 
for n-octane except that the air pressure was kept constant at 1 atmos- 
phere and the air temperature varied. Figure 9(c) shows temperature- 
time histories for n-octane where all the droplets had a 50 -iJiicron 
initial radius and both the air temperature and pressure were varied. 

In this case none of the nine temperature -time histories shown would 
have the same asyrr^totic temperature. The end lines joining the curves 
at the same air temperature represent, therefore, a variation of asymp- 
totic temperatures with air pressures at three different air temperatures. 

Figures 10(a) to 10(c) and ll(a) to ll(c) show the same histories 
described in figures 9(a) to 9(c) but for benzene and n-pentane^ 
respectively. 

Figure 12(a) shows temperature-time histories of the three fuels 
where the initial droplet radius was 50 microns, the air pressure was 
1 atmosphere, and the air temperature was varied. Figure 12(b) shows 
the three fuel droplets having the same initial radius, with the air 
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■bcmperat/Ure kept const-ant at 500^ F and the air pressure varied. It 
am be seen from these curves tliat the lover volatility fuels have 
higher asymptotic temperatures and take a longer time to reach them for 
the same operating and ambient conditions. By plotting the results of 
the present calculations^ it vas found that^ with ail other conditions 
constant^ the time for droplets to reach the vet-bulb temperature varied 
approximately as the initial radius to the 1.8 power. 

In the combustion chambers the mass of the fuel droplets vaporized 
and the penetrations of these droplets relative to the air are essential 
in deciding the zone where combustion is likely to take place. Fig- 
ures 15(a) throiigh 15(i) show time histories of droplet penetration d, 
relative to the air, and mass, percent of original, M during the 
unsteady state for n-octane at nine different ambient air conditions . 

Each of these figures contains the previous information for the three 
initial droplet sizes and two initial droplet velocities. As the tem- 
perature decreases but with the same air pressure, the droplets take 
a longer time and penetrate farther before reaching their steady states j 
as the air pressure increases but for the same air temperature, the ^ 
droplets take a longer time but penetrate less before reaching their 
steady states. 

Figure Ik presents a comparison of benzene and n-octane vet-bulb 
temperatures calculated in this report with vet— bulb temperatures 
experimentally obtained by Ingebo (ref. 2). Referring to the benzene 
curves, it can be seen that the calculated vet-bulb temperatures are 
always somewhat higher than the experimental vet-bulb temperatiores . 

This would indicate that the computed mass transfer vas too low or the 
computed heat transfer was too high or both. Referring to the n-octane 
curves, it can be seen that there is a cross-over of the two curves. 
Extrapolation of the data of Ingebo would indicate that the droplet 
would boil, in contradiction to the diffusion theory. Whether this is 
a deficiency in the theory or in the experimental data of Ingebo is not 
clear. 

A typical temperature -time history obtained in the experimental 
investigations is shown in figure I5 for a droplet of n-octane of approxi- 
mately 750-micron initial radius when exposed to air moving at a constant 
velocity of about 16 fps. The portion 1-2 of this graph shows the output 
of the thermocouple when initially in room air. At point 2 the droplet 
was hung on the thermocouple. Portion 2-5-4 of the graph represents the 
rapid response of the thermocouple to the cooler-than-room-air droplet 
followed by a gradual cooling to the wet -bulb temperature corresponding 
to room air at point 4. This period is of particular interest in that 
it gives an indication of the speed of response of the thermocouple and 
recording instrument. 

At point 4 a stream of hot air was suddenly swung under the droplet. 
The unsteady-state or heating-up portion of the curve can be seen from 
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point 4 to a point chosen arbitrarily at 5 where the droplet would be 
judged to have closely reached the wet-bulb terrperature corresponding to 
the temperature of the air stream. The portion 5-6 of the graph is con- 
sidered to be the steady-state vaporization period. The last part of 
this portion, near 6 , is of questionable reliability, since observation 
showed that as the droplet nears complete vaporization it deviates from 
the spherical shape assimed in the calculations and, instead, takes a 
distorted shape depending on the exact shape of the thermocouple junc- 
tion. It should also be noted that the purity of the fuels was only 
about 98 percent and that this portion of the curve may represent the 
vaporization of the high-boiling -point impurities in the mixture. 

At point 6 a thin film of the liquid fuel still remains on the 
thermocouple element. This soon disappears and the thermocouple tem- 
perature rises to the air temperature as represented by point 7. At 
that point the thermocouple is completely dry and the temperature 
recorded at 7 is the air temperature dioring the entire plot. Since 
this takes only a short time, the air temperature is not likely to have 
varied during this period. This procedure has the advantage of using 
a single thermocouple to measure both the droplet and air temperatures . 

Since the amount of liquid fuel remaining as a thin film on the 
thermocouple element at point 6 in the graph is minute compared with 
the total liquid contained in the initial droplet, the relative magni- 
tudes of time of portions 4-5 nnd 5"6 of the graph can be assumed to 
represent the relative importance of the unsteady and steady states, 
respectively. 

Figure 16 shows several time histories of droplets of n-octane all 
of approximately equal size as subjected to air of varying temperature. 
The portions of the curve corresponding to 2-5-4 in figure I5 have been 
omitted here. It can be seen from these curves that the unsteady-state 
portion of the vaporization time assumes greater importance and becomes 
dominant at higher air temperatures. 

Figure I7 shows time histories of liquids of varying volatility 
vaporizing under the same conditions of pressure, temperature, and 
velocity. 

gpjice in the calculated curves presented in this report the velocity 
oP the droplet was made to decrease relative to the air because of aero- 
dynamic drag forces, while the experimental curves were measured at 
essentially constant air velocity past the droplets, caution must be 
exercised in comparing the experimental results thus obtained with the 
calculated curves. Figure I8 shows two calculated temperature, velocity, 
and rate of mass transfer (unsteady -state) histories of two droplets of 
n-octaae of 50-micron initial radius and 100-fps initial velocity each. 
One was made to slow down, while the other remained at 100-fps velocity 
throughout the calculated period. Other conditions were kept equal. 
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It can be seen that the constant-velocity droplet spends less time in 
the unsteady-state period but has a higher rate of mss transfer during 
that period with consequently shorter total vaporization time. The n^ 
result would be that the unsteady state is relatively longer for the 
cLToplet that is moving at a constant speed. 


The following deductions can now he made. 


(1) For the 
in magnitude with 
air temperature . 


same fuel the unsteady state becomes relatively larger 
respect to the total vaporization time the higher the 


(2) For the same air temperatures the unsteady state is relatively 

vapo^zatlon tl« the higher 

( 3 ) High- volatility fuels have lower wet-bulb ten5)eratures for the 
same air conditions. 


( 4 ) At relatively low air temperatures extremely low volatility 
fuels have wet-bulb temperatures close to the air temperatures and spend 
only a vepr small portion of their vaporization time in the unsteady 
state. At high air temperatures, however, the difference between the 
air and wet-bulb temperatures increases. 


(5) For any one fuel the wet-bulb temperature is higher the hicher 
the air pressure or temperature. 

penetration - mass -vaporized relationships of different size 
droplets ^e of particular interest in jet-engine combustion-chamber 
esign. When the fuel is injected continuously into a stream of moving 

injector nozzle first as a ligament or sheet after 
which this ligament or sheet breaks down to different size droplets 
originally movi^ at the same speed. If the velocity of these droplets 
s higher than the air velocity, these droplets start slowing down rela- 
tive to the air, with the smaller droplets slowing down faster, and 
penetrate a lesser distance while losing their mass faster and earlier 
a ong their travel path. The mass of vapor given away by each droplet 
IS carried along with the air through the combustion chamber. To this 
continuously added the vapor given away by the larger droplets 
ttet have been formed later in time but that reached that point beLuse 
of their more rapid and greater penetration. The total amount of fuel 
vapor present at a certain cross section of the combustion chamber there- 
fore consists of the total vapor given away by all the drops before 
reaching tlmt cross section. In order to have a short combustion cham- 
ber, a combustible air-vapor mixture should be formed a short distance 
after Injection. This would necessitate that the largest number of the 
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drop sizes fall in the smaller radius range so that the weight distri- 
bution of different droplet sizes in the spray would not he weighted too 
highly in favor of the droplets of large radius range, since these large 
droplets have a much greater penetration before they contribute much 
vapor toward forming a combustible mixture. 

To illustrate the possible use of the ciorves of droplet mass versus 
time for combustion-chamber design purposes, and to empliasize the impor- 
tance of the unsteady portion of the time history of the droplets, the 
following example was worked out with certain arbitrary but reasonable 
assumptions made as to the combustion— chamber operating conditions . 

These assumptions were: 

(1) No effect of one evaporating droplet upon the other. This is 
in obvious contradiction to the facts, but since this effect is unknown 
it was necessary to neglect it. 

(2) A constant air velocity of 5O fps. 

(5) An air pressure of 1 atmosphere and a temperature of 500° F. 

The 500° F air temperature is admittedly high with respect to the air 
temperature coming from the compressor, but calculations were available 
for this temperature and recirculation, which was neglected, would 
undoubtedly increase the air temperature. 

(1+) The fuel used was n-octane. 

(5) A drop-size weight distribution of 25 percent for the 10-micron 
group, 65 percent for the 50-micron group, and 10 percent for the 
100-micron group. 

(6) An initial fuel velocity relative to the air of 100 fps. 

The previously performed calculations were extended through the 
steady state until the droplets nearly disappeared.' Such a calculation 
showing temperature, mass, and penetration relative to the air versus 
time for the 50-micron droplet is shown in figure 19 . The masses vapor- 
ized from each size droplet were then calculated and plotted versus dis- 
tance traveled by each droplet relative to the combustion chamber using 
the assumed airspeed of 50 fps, as shown in figure 20. 

Using the arbitrarily chosen droplet size distribution, the percent- 
age of the total fuel vaporized from each droplet size group up to any 
distance after droplet formation was then plotted in figure 21. The 
total vapor present, which is the sum of the three curves, is also shown 
in figure 21. Assuming a fuel-air mixture strength in the primary com- 
bustion zone to be in the neighborhood of O.O8, and a possible inflamma- 
bility ratio of O.O5 of the vapor-air ratio, a total vaporization in the 
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neighborhood of 35 percent would be adequate for initiating combustion 
^der these assumptions. According to the calculations presented in 

figure 21, this would take place 8 inches after the droplets have been 
I orinea . 


It IS to be noted that these calculations are approximate at best, 
■pjrbulence, backflow, and heat transfer by radiation have not been taken 
into account. However, they serve to show the possible use of the cal- 
culated ci^es for estimation purposes. A knowledge of droplet size 
distribution seems also to be fundamental in any such vaporization 
calculations . 

An important observation shown by the above calculations is that 
^ e 5 O- and 100-micron droplets reach the combustion zone while still 
in the heating -up or unsteady-state portion of vaporization, a fact 
that further emphasizes the importance of the unsteady state in anv 
combustion studies. It is the opinion of the authors that the unsteady 
state IS of primary importance in forming a combustible mixture and that 
experimental studies which could be used either to verify or to correct 
the extrapolation of the heat and mass transfer correlations used should 
e conducted. Until these extrapolations are confirmed, the absolute 
acc^acy of the curves can be questioned, although it is not believed 
ttet any of the conclusions reached would be altered. After these rela- 
tionships have been confirmed or corrected, it would then seem necessary 
to determine the effect of one droplet upon another. This information 
together with drop size distribution should then permit a rational and 
reasonable determination of the time required to form a combustible mix- 
ture in a heterogeneous mixture of fuel droplets and air. 


University of Wisconsin, 

Madison, Wise., March 12, 1953. 
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APPENDIX A 

SAT^LE OF CALCULATIONS 


The unsteady-state portion of the vaporization time of a droplet of 
benzene was calculated for the following data: 

Droplet initial radius - 50 microns = 1,9685 X 10“5 in. 

Droplet initial temperature = 50^ D = 510^ R 

Droplet initial velocity relative to the air = 400 fps = 

4,800 in. /sec 

Air pressure P^i = 1 atm = 14.696 Ib/sq in. abs 

Air temperature Tg = 1,000^ F = l,46o° R 

A stepwise method was used in integrating the equations as follows: 

The temperature of the droplet was assumed to undergo a small increase ATg 
over a small period of time A0, during which the properties of the liquid 
fuel, fuel vapor, and air could be considered to be constant. The heat 
and mass transfer during this period were then calculated together with 
the velocity change, from which the time required for the temperature 
increase together with the radius and velocity at the end of the increment 
were calculated. However, these three values had to be assumed prior to 
the calculation and compared with the calculated values . 

For a droplet -temperature increase of 5^ ^ assumed that: 

= 2.7 X 10-5 sec 

^2 = ^>725 in. /sec 
Tq 2 = 1-970 X 10"5 in. 


For this increment: 


= average temperature of droplet = 


510 + 515 


= 512.5 R 


2 
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The various properties in the film and for the droplet were then 
calculated and read as follows: 

p^ = 0.970 Ib/sq in. abs 
Pjjj = 2.678 X 10“5 Ib/cu in. 

Pm = 1.475 X 10-6 lb/(in.)(sec) 

= 5-591 X 10-7 Btu/(in.)(sec)(°F) 

Cp^ = 0.2609 Btu/(lb)(°F) 

- 3-60 X 10 “^ sq in. /sec 
A = 189.0 Btu/lb 
Cpp = 0.422 Btu/(lb)(^) 
ppp = 3.15^ X 10 “^ Ib/cu in. 

Pj- 2 = 3 -150 X 10“^ Ib/cu in. 

Cp^ = 0.4080 Btu/(lb)(°F) 


1*0 = Average radius of the droplet for the increment 

= ^ol ^02 
2 

= ^ 1.970j ^ ,0-5 

= 1.9693 X 10"^ in. 
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U ^ Average velocity of the droplet 

_ Ul U2 
2 

- + 4,72^ 

2 

= 4,762.5 in, /sec 

Nj^e = Reynolds number 


= 540.568 

Cd (^6 Npjgj 
Applying equation (18): 

^ = _lc ^ 

^ ® ^ 9li ^o 

and 

U = 65.681 in. /sec 


= 0.6633 


2.4526 in./sec^ 


Therefore 


U2 = 4,800 - 65.681 = 4,734.319 in. /sec 


as compared with the assumed value of 4^725 in. /sec. This process was 
then repeated^ selecting the correct coefficient of drag each time until 
^ for was arrived at that satisfied the velocity eq^uation. 

Therefore 


^2 ^ ^j» 730.79 in. /sec 


32 


NACA TN 3179 


and 


%e = 3^^0.775 


Therefore 


N, 


Sc 




= 1.5299 


kn- = 


D p 
v'^m 


2Pa^o 


2 + 


= 2.51^27 X 10 “^ sec -1 


a = 


PfL 


= 1.050 


logp 


~ PfL 


Therefore 


w = kgapfLUitr^^ = 1,2381 x 10"^ Ib/sec 


Therefore 


Cpm^^m 


h = 


hL 

2r, 


o *— 


Npj. = = 0.6883 


2 + 0-6(%r)"''^(l»Re)'''^ 


= 1.6721 X 10“5 Btu/(sec)(sq in.)(°F) 


Q = hAo(Tp - Tl) = 7.7213 X 10-5 Btu/sec 


- 1 


= 0.9683 
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Therefore 

Q = Q E = 7.4T6i<- X 10~5 Btu/sec 

^ - 1 

Qg = Q - Qy. = 2A488 X 10"^ Btu/sec 
= wA = 2.54006 X 10"^ Btu/sec 

and 

Ql = Qy - = 7.2424 X 10"5 Btu/sec 

= Mass of droplet at beginning of increment 

= 510° R) 

= 1.00956 X 10“9 lb 

m 2 = Mass of droplet at end of increment 
= m^_ - w A0 
= 1.00901 X 10"9 lb 

m = = 1.009186 X 10"^ lb 

2 

Therefore 

mCpT ATl 

A0 = = 2.8426 X 10"^ sec 

Ql 


(as coinpared with the estimated value of 2.7 X 10“^ second). 


NACA TN 5179 


3 ^ 


I 3"l2 , 

“■02 - "■ 

(as con 5 )ared with the estimated value of 1.970 X 10"5 inch). 

The calculations for this first Increment were then repeated taking 
the calculated values of time, velocity, and radius as new assumptions 
to obtain more accurate results. The final values of U2, 

calculated for the first increment, then beceune the original values of 
the next increment and the same stepwise method was used for the second 
increment. The calculations were stopped when Ql approached zero, 

indicating that the asymptotic or wet-bulb temperature has been approached. 
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APPENDIX B 

CALCULATIONS OF PROPERTIES IN FILM 


The comp'u'ta'bions described in the analysis req^uire the use of the 
thermal conductivity, viscosity, density, and specific heat of the air- 
vapor mixture in the film surrounding the droplet. The concentrations 
of fuel vapor and air are not constant across the film, with the mixture 
being the richest in the vicinity of the droplet and leaner away from 
it. To a good approximation the modal ratio of fuel vapor to air equals 
the ratio of the partial pressure due to the fuel vapor to the partial 
pressure due to the air in the film. Thus the largest fuel-vapor - air 

ratio occurs at the surface of the droplet and is essentially equal to 

the ratio of the vapor pressure at the droplet ten5)erature to the partial 

pressure due to the air, or Pfp/PaL* This ratio becomes zero at the 

outer edge of the film, provided there is no interaction between droplets. 

It has been shown for plane surfaces (ref. 8 ) that in the case of a 
fully permeable film, where mass transfer takes place in both directions 
and equimolal diffusion exists, the partial pressure gradients would be 
linear throughout the film. The temperature in the film which in most 
of the cases varies from a maximimi of Tp, the air temperature at the 

outer edge of the film, to Tp, the droplet temperature at the inner 
edge, would follow a linear gradient only in the case of no mass trans- 
fer. In the case under consideration, that is, a semipermeable film 
with mass transfer in one direction only, where part of the incoming 
heat is picked up by the diffusing vapor, the linearity of both the ^ 
partial pressure and temperature gradients in the film ceases to exist. 

An exact evaluation of the properties in the film surrounding a 
spherical droplet would render the equations almost impossible to inte- 
grate, and furthermore the precision with which some of the properties^ 
are known probably does not justify the work required. Thus an averaging 
technique in the film was employed. The properties then have been cal- 
culated at the logarithmic average temperature in the film T^ and for 
a fuel-vapor - air mixture of constant concentration equal to one-half 
the concentration at the droplet surface. The following equations were 
used: 


T, 


m 



loge 


Tb 


(29) 
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(50) 



^fL\ PfL 


12RT^ 



(31) 


(52) 



In these relations the pure properties of the fuel vapor and air 
have been taken at the logarithmic mean temperature in the film. 

Figure 22 shows the properties in the film versus droplet tempera- 
ture as calculated for n-octane and air. 
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APPENDIX C 

JUSTIFICATION FOR ASSUMPTION OF HIGH FUEL CONDUCTIVITY 
AND UNIFORM DROPLET TEMPERATURE 


As can be deduced from equation (2l), the calculations performed 
In this project assumed a large thermal conductivity through the liquid 
droplet and a resulting uniform temperature throiaghout the droplet at 
all times. An estimation of the validity of the above assumptions can 
be made from the charts prepared by Gurney and Liirle (ref. 9j discussion 
given in ref. lO) for determining the relationship between temperature 
and time at various points in spheres and other shapes . To adapt this 
method to the cases of interest here, the resistance ratio ^l/^L^o 

used, where is the thermal conductivity of the liquid fuel evaluated 

at the temperature of the droplet; hL. a coefficient of heat transfer to 

the liquid droplet based on the portion of heat transferred to the droplet 
that appears as sensible heat (Ql)j radius of the droplet. 

The use of this parameter together with the charts permits an estimation 
of the possible difference between the temperature at the surface of the 
droplet and the temperature at its center^ respectively. The Gurney- 
Lurie charts specify that the temperattires at the midpoint and at the 
surface would differ considerably when the resistance ratio is small and 
only slightly when it is large. Referring specifically to the Gurney- 
Lurle charts for spheres, if the resistance ratio is 6 or more, the tem- 
perature difference is negligibly small; while, if the resistance ratio 
is 1 or 2, the temperature difference is appreciable and should be taken 
into account. However, the Gurney-Lurie charts were based on the assump- 
tions of constancy in the temperature surrounding the surface of the body, 
the resistance ratio, and the thermal diffusivity in the body under con- 
sideration; and it should be realized that, under the conditions of an 
evaporating droplet, the charts should be used for estimation piorposes 
only. 


Values for the resistance ratio are plotted in figure 25 Tor n-octane, 
benzene, and n-pentane versus heating -up time together with the droplet 
temperatures calculated by the methods previously described in this report 
for an air pressure of 1 atmosphere and an air temperature of ^ 00 ^ F. 
Figures 24(a) and 24(b) show the resistance ratio plotted for benzene at 
a constant air pressure of 1 atmosphere and a constant air teirperature 
of 500 ^ F, respectively. 

The data of figures 25 and 24 indicate that the assumption of high 
thermal conductivity in the liquid droplet is unquestionably good for 
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liigh-volatilitiy fuels ^ is of reasonable validity for mediuin-vola’fcili'by 
fuels ^ and should be further investigated for the low-volatility fuels. 
If internal circulation (ref. 7) exists in t' : droplets under the con- 
ditions of interest here^, it could effectively increase the conductivity 
of the liquid. 
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TABLE I 

REFERENCES FROM WHICH PROPERTIES WERE OBTAINED 
FOR FUELS AND FOR AIR 


Property 

Reference 

n-Octane 

Benzene 

n-Pentane 

Air 

Pl 

11 and 12 

11 

13 

— 

kf 

11 

12 

11 and 12 

1 

Pf 

11 

12 

11 

1 

A 

11 

12 

11 

— 

PfL 

11 

11 

11 

— 

Cpf 

14 

i4 

i4 

11 

^pl 

13 and l4 

15 

13 and l4 

— 

^L 

15 

15 

15 

— 
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Figure 2.- Typical vaporization-time history for droplet of n-octane 
with initial radius of 50 microns. T^_ = 50 ^ = ^00 fps; 

Prp z= 0.5 atmosphere; T-d = 1,000° F. 
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(a) Effect of droplet velocity, = 10 microns. 



(b) Effect of droplet radius. Velocity, 100 fps. 

Figure 5.- Effect of droplet velocity and radius on T„^. 

= 1 atmosphere; = 100° F. 
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Figure 10,- Continued. 
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Figure 9*- Continued 
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(c) With constant initial radius of 50 microns 







= 500 



Figure 11.- Concluded. 
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Figure 12.- Concluded. 
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(a) Pip = 0.5 atmosphere; Tg = 100° F. 

Figure 13.- Penetration-time and mass-time histories during vinsteady 
state for n-octane at nine different ambient air conditions. 
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(b) 


Prp = 0.5 atmosphere; Tg = 500° F. 


Figiire 13.- Continued. 
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(e) Fiji = 1 atmosphere; Tg = 


Figure 15 •- Continued 
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atmospheres ; Tg = 500° F. 


15. - Continued 
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(i) Prj, = 5 atmospheres; Tg = 


Figure 15 •- Concluded 










Temperature 



Figure I 5 .- Temperature-time history of a droplet of n-octane. 

1- 2 Room air temperature 

2 Droplet placed on thermocouple 

2- 3 Rapid self-cooling to approximately wet-bulb temperature 

at room air temperature 

3 - k Gradual approach to wet-bulb temperature 
4 Hot air suddenly swung under droplet 

4- 5 Heating-up or unsteady-state period to wet-bulb temperature 

corresponding to hot air 

5 - 6 steady state at wet-bulb teniperat\ure 

6- 7 Heating up of thermocouple to air temperature 



Temperature 


74 


NACA TN 3179 



Figure I6.- Effect of air temperature on heating up of n-octane. 


Temperature 
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Figure 17*- Effect of fuel conipositlon on heating-up time at constant 

air temperature. 
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Figure 18.- Tenrperature^ velocity, and rate of masc transfer histories 
of two droplets of n-octane of 50-micron initial radius. Tj^ = 50° 

Uf = 100 fps; = 1 atmosphere; Tg = 1,000*'* F. 



Figure I9.- Temperature- mass-^ and penetration-time histories of 

50-micron n-octane droplet. 
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Figure 21.- Mass vaporized of each droplet size group and total mass 
vaporized (percentage of total fuel injected) . 




Droplet temperature, ®F 


Figure 22.- Average properties of air- vapor mixtiore in film at 
different droplet temperatures. Air pressure, 1 atmosphere 
air temperature, 200° F. 
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(a) At constant air pressure of 1 atmosphere. 

J^igure 24.- Resistance ratio and droplet temperature for henzene. 
Up = 100 fps; = 50 microns. 



Resistance ratio, k^/h^rQ , unitless 
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Time, milliseconds 

(b) At constant air temperature of 500° F. 
Figure 2k.- Concluded. 
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